Electromagnetic Energy, Absorption, and Casimir Forces. I. Uniform Dielectric Media 

in Thermal Equilibrium 
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The derivation of Casimir forces between dielectrics can be simplified by ignoring absorption, cal- 
culating energy changes due to displacements of the dielectrics, and only then admitting absorption 
by allowing permittivities to be complex. As a first step towards a better understanding of this situ- 
ation we consider in this paper the model of a dielectric as a collection of oscillators, each of which is 
coupled to a reservoir giving rise to damping and Langevin forces on the oscillators and a noise po- 
larization acting as a source of a fluctuating electromagnetic field in the dielectric. The model leads 
naturally to expressions for the quantized electric and magnetic fields that are consistent with those 
obtained in approaches that diagonalize the coupled system of oscillators for the dielectric medium, 
the reservoir, and the electromagnetic field. It also results in a fluctuation-dissipation relation be- 
tween the noise polarization and the imaginary part of the permittivity; comparison with the Rytov 
fluctuation-dissipation relation employed in the well-known Lifshitz theory for the van der Waals 
(or Casimir) force shows that the Lifshitz theory is actually a classical stochastic electrodynamical 
theory. The approximate classical expression for the energy density in a band of frequencies at 
which absorption in a dielectric is negligible is shown to be exact as a spectral thermal equilibrium 
expectation value in the quantum-electrodynamical theory. Our main result is the derivation of an 
expression for the QED energy density of a uniform dispersive, absorbing media in thermal equilib- 
rium. The spectral density of the energy is found to have the same form with or without absorption. 
We also show how the fluctuation-dissipation theorem ensures a detailed balance of energy exchange 
between the (absorbing) medium, the reservoir and the EM field in thermal equilibrium. 

PACS numbers: 42.50.-p, 03.70.+k, 



I. INTRODUCTION 

Based on the assumption of an electromagnetic (EM) 
energy \hu) per mode of angular frequency u> at zero tem- 
perature, Casimir [l[ famously showed that there is an 
attractive force between two uncharged, perfectly con- 
ducting plates. Lifshitz generalized the theory to the 
case of two thick dielectric slabs in thermal equilibrium 
by calculating the stress tensor for the fluctuating field in 
a vacuum between the slabs. Casimir's original method 
involving changes in zero-point field energy was later ex- 
tended to dielectrics by van Kampen et al. Q in the 
"nonretarded" case of small separations, and by others 
for arbitrary separations 0]. 

Derivations of the Lifshitz formula that invoke changes 
in zero-point energy begin by assuming real dielectric 
permittivities. After an integral over frequencies for the 
force function of the distance separating the di- 

electrics is obtained, the permittivities are allowed to be 
complex functions, analytic in the upper half of the com- 
plex frequency plane as required by causality. This allows 
an analytic continuation to an integral involving only 
purely imaginary frequencies, at which the permittivities 
are purely real, and the resulting expression is equivalent 
to that of Lifshitz, who requires complex permittivities 
(absorption) through the fluctuation-dissipation relation. 
In this paper we take a first step towards a better under- 
standing of why this approach leads in the end to results 
that are equivalent to those obtained by Lifshitz-type ap- 
proaches, which are not based on calculations of energy 
changes and which account explicitly for absorption via 



a fluctuation-dissipation relation. 

We begin in the following section by revisiting the clas- 
sical theory of electromagnetic energy density for quasi- 
monochromatic fields in a dispersive, absorbing medium. 
At frequencies u> at which absorption is negligible the 
classical expression for the average energy density is @ 
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where e R and [ir are respectively the (real) permittiv- 
ity and magnetic permeability at frequency u), E w (r) 
and H w (r) are the electric and magnetic fields at u> and 
the bar over u(r, ui) indicates a time average (see sec- 
tion II. A). In the following section it is also shown that, 
within a band of frequencies for which absorption is neg- 
ligible, u(r, u) gives exactly the spectral energy density as 
long as the fields at different frequencies undergo uncor- 
rected fluctuations. In the quantum-electrodynamical 
(QED) theory in which there are quantum field fluctua- 
tions, this holds for the expectation value of the thermal 
equilibrium energy, and in particular for the zero-point 
energy. More interesting for our purposes, however, is 
the fact that these expectation values in QED have the 
same form with or without absorption in uniform media. 

This paper is organized as follows. In Section [IT] we 
review some aspects of the classical theory of electro- 
magnetic energy density in dispersive, absorbing media. 
Equation UJ is shown to give the total energy density 
of the field and the polarizable particles of a purely di- 
electric medium (/i = 1) at frequencies for which dissi- 
pation is negligible. Whereas it gives approximately the 
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total energy density in the case of quasi-monochromatic 
fields, it gives exactly the average total spectral energy 
density within a band of frequencies in which different 
frequencies undergo uncorrelated fluctuations, provided 
that absorption can be ignored in this band. We review 
in Section IHII the total gMcmiitm-electrodynamical ther- 
mal equilibrium energy density for a dispersive magnc- 
todielectric medium in which absorption is negligible. In 
Section HVl we treat in some detail the quantum theory of 
a dielectric medium modeled as a collection of polarizable 
material oscillators, allowing for dissipation by coupling 
each of these oscillators to a reservoir and deriving the 
fluctuation-dissipation relation for the noise polarization 
arising from quantum fluctuations of the bath oscillators. 
We argue, based on this derivation, that the Lifshitz the- 
ory employing such a fluctuation-dissipation relation is in 
fact a classical stochastic clectrodynamical theory. Pro- 
ceeding as in the Lifshitz theory in which the noise polar- 
ization acts as a source of a fluctuating electromagnetic 
field, we obtain expressions for the quantized electromag- 
netic field in a dissipativc medium; these quantized fields 
have the same form as in the Huttncr-Barnett theory in 
which the Hamiltonian for the coupled system of oscilla- 
tors for the dielectric medium, the reservoir, and the elec- 
tromagnetic field is diagonalized [(| . Adding each contri- 
bution to the total energy density, including that from 
the reservoir, we derive (for the first time to our knowl- 
edge), the expression (|106p for the QED energy density of 
a uniform dispersive, absorbing medium in thermal equi- 
librium. We show in particular that it has exactly the 
same form as obtained in Section [IIII when dissipation is 
neglected, and is affected by dissipation only through the 
dependence on dissipation of the real part of the refrac- 
tive index. In Section [V] we summarize our conclusions 
and discuss briefly how they apply to two examples: the 
Einstein A coefficient for spontaneous emission and the 
van dcr Waals interaction of atoms embedded in a dissi- 
pativc dielectric medium. 

This is intended to be the first of two papers dealing 
with electromagnetic energy in dispersive, dissipative me- 
dia. The approach described here will be extended in a 
forthcoming paper on effects of dissipation when Casimir 
effects (presence of material boundaries) are calculated 
following an approach based on zero-point energy [?J . 



II. CLASSICAL ELECTROMAGNETIC 
ENERGY DENSITY IN DISPERSIVE, 
ABSORBING MEDIA 

The problem of defining and calculating an electro- 
magnetic energy density for the classical EM field in dis- 
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persive and absorbing media was investigated in detail 
by Barash and Ginzburg [H, Q . Here we review some of 
the main points and difficulties associated with this prob- 
lem. The classical expression for electromagnetic energy 
density can be derived from the Poynting theorem, 
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where S = c (E x H) /An is the Poynting vector in the 
conventional notation. We assume that the constitutive 
relations connecting E, D, B and H are linear, isotropic 
and spatially local, so that, writing 
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We write the (complex) electric permittivity e and mag- 
netic permeability \i in terms of real and imaginary parts: 
e(r,w) = e R (r,u) + iej(r,u;) and (j,(t,ui) = Hr(t,u)) + 
«/x/(r, tt>). At this point we should stress that we con- 
sider only passive media throughout this work, meaning 
that e/(w), /i/(o;) > for all frequencies. The well-known 
relations 



E(r,-w)=E*(r,w) , B(r, -u) = B*(r, u), 
e(r,-w) = e*(r,w) , fi(r, -w) = fi* (r, u) (5) 

follow from the reality of E and B. Then 
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obtained using the properties §5§. (To simplify notation 
henceforth we do not indicate any dependence of e and 
/i on r.) The constants of integration vanish under the 
assumption that E(r, t) — » and H(r, t) — > as i — ► — oo. 

So far we have been using the concept of electromag- 
netic energy density a bit loosely, and at this point we 
would like to make our statements more precise. Being a 
direct consequence of Maxwell's equations, the balance 
relation ^ is valid under arbitrary thermodynamical 
conditions, so it can be used to describe general out- 
of-equilibrium systems. However, in those situations one 
should be careful on defining energies, as it is easy to 
see that in such systems the r.h.s. of (J5]) contains not 
only the "standard" electromagnetic energy W(r, t) but 
also the dissipated heat Q(r,t), and, as it is discussed 
in length by Barash and Ginzburg 0, in general it 
is impossible to separate the two in a unambiguous way. 
So, in order to avoid confusion, we always work with the 
sum 



W{v,t)=W{v,t)+Q{v,t), 
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and we refer to W(r, t) as the electromagnetic energy for 
the lack of a better term, but always bearing the above 
considerations in mind. In the case where thermal equi- 
librium is established, the evolved heat Q(r,t) vanishes 
and W(r, t) coincides with W(r,t). In addition, even 
for situations in thermal equilibrium, we shall make the 
distinction of the EM energy when there is absorption 
present (e/,^/ ^ 0), that we will call W(r,t), and when 
absorption is absent (ej = fii = 0), that we will denote 
by u(r, t). 

In the limiting case of a monochromatic field we have 



E(r,w) = -E Uo (t)[6(lo - lo ) + 5(lo + lo )], (10) 



and, if absorption at frequency loq is negligible, it follows 
from ([7]) that at thermal equilibrium 



A. Quasi-monochromatic fields 

In the case of quasi-monochromatic fields it is conve- 
nient to write the fields E(r, t) and H(r, t) as 

E(r,t) = i [EofotJe-^+ESCr.t)^ 1 "*] 

= / c^[E(r,u;)e- <wt + E*(r,u;)e iwt ], (13) 
Jo 

H(r,t) = i [HoCr^e-^+HSCr,^*] 

= / duj[H(r,uj)e~ wt +H*{r,w)e iut ], 
Jo 

(14) 

where the envelope functions Eo(r, t) and Ho(r, t) vary 
slowly in time compared to e~ lulat . The Fourier compo- 
nents E(r, lo) and H(r, to) in this case are sharply peaked 
at the frequency loq, and we consider the time-averages 
{We and Wh) of We and Wh over times long compared 
to 2tt/loq but short compared to times over which E and 
Ho vary significantly. We also assume that e(tu) and /Lt(cj) 
vary slowly enough near w = wo so that we can retain 
only the first-order terms in their Taylor series about loq 
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With these approximations we find straightforwardly 
that 
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where we have averaged over the period 2tt/loq and have 
used 
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for e(jjo) = cr(lo) and lo 1 — > lo. Together with the corre- 
sponding result for Um{ t ), this gives Eq. |T|). 
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where we have used 
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The time t here has been assumed to be short compared 
to the time over which the slowly varying field envelopes 
Eo(r, t) and Ho(r, i) change significantly, as otherwise 
even small deviations from monochromaticity can inval- 
idate the approximations (JT6J) and (fT7|) . 

Equation (|16|) gives the first two terms corresponding 
to the expression (8) of Barash and Ginzburg Q. To 
obtain the remaining terms in that expression we must 
include terms proportional to dej/du, that may be writ- 
ten as 
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and, after straightforward manipulations, as 
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which is nothing but the time integral of the third 
term in Eq. (8) of Q. Since the discussion of quasi- 
monochromatic fields is not the main goal of this paper 
we shall stop it here, referring the reader to 0, 0] for 
further details. There one can find extensive discussions 
on the interplay of EM energy and heat generated, on 
how dissipation makes it (in general) impossible to write 
the "standard" EM energy W(r, t) in terms of dielec- 
tric functions alone, how it is possible to go beyond the 
quasi-monochromatic approximation when dissipation is 
absent, etc. 



B. Uncor related frequencies 

Another situation where it is possible to simplify the 
general expression for the energy density is when we have 
stochastic fields such that their autocorrelation functions 
in the frequency domain satisfy 

«E(r, u) • E*(r, u'))) = \ |E(r, u)\ 2 5(u' - u), 

«H(r, u) • H*(r, u'))) = i|H(r, u)\ 2 6(u> - u), 
((E(r,w) -E(r, w ')» = «H(r,w) • H(r,a/))) = 0, (21) 



where ((...)) denotes the average over the appropriate 
ensemble. From these correlation functions it follows that 
correlations in the time domain are stationary: 



«E(r,i)-E(r,i')» = F E (v,t-t'), 
«H(r,i).H(r,i')» = F H (r,t-t'). 



(22) 



It is clear that equations ([22)) arc satisfied when there 
is no net dissipation or gain, which means that cither 
ej = jjbi = or that there are Langevin-type forces in the 
system that compensate for dissipated energy. Restrict- 
ing ourselves to the first case in this simple example, we 
can use (|2"Tj) to calculate the ensemble average of ([7]) and 
(jSJ) and obtain at thermal equilibrium 
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Of course this expression is strictly valid only over fre- 
quency ranges at which absorption is negligible; in such 
ranges the integrand of Eq. ([23]) gives exactly the spectral 
energy density. 



C. Classical oscillator model for the energy density 

In order to better focus on the physics involved in these 
considerations of energy density we briefly review a clas- 
sical model in which the medium consists of N har- 
monic oscillators per unit volume, each having a natural 
oscillation frequency Q and satisfying the equation of mo- 
tion 



x + fTx = — E. 

m 



(24) 



(For notational simplicity we do not indicate here the r 
dependence of E.) The polarization density and dielectric 
constant in this model are respectively 



P = iVex = 
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where u p = (AirNe 2 /m) 1 / 2 is the plasma frequency. We 
write Poynting's theorem in its integral form: 
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The integral of the normal component of S on the left- 
hand side is, as usual, over a surface enclosing a volume 
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V. From (HI) 



is 
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and therefore 
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u is the density of total energy, that in the field plus that 
in the medium. Using 

elm . uje/m 

x = 7^ ^E^coscjt and x=- — - -E^smwi, 

\V — UJ 1 \V — uj z 

(30) 

and (|26| for a monochromatic field E^coswi, we find 
after cycle-averaging that 
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confirming that equation ((T|) defines the totaJ energy den- 
sity of a dielectric medium (fin = 1). From the relation 
H 2 = ei?(a>)E 2 , 
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the term (l/167r)w[de i? /dw]E 2 is seen from ([29]) and ([30]) 
to be the (cycle-averaged) kinetic energy per unit volume 
of the material oscillators in this model. 

Absorption is included in this model by adding Tx (T > 
0) to the left-hand side of flSJ) [13 ■ Then it is easily 
shown that the rate of change of energy density W in the 
volume V, defined such that 
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A similar expression is derived in the QED theory in 
Section Hvl An important difference, however, is that in 
the QED theory there is an additional term arising from 
Langevin forces, which are required for the preservation 
of commutation relations. In addition, we should stress 
that these Langevin forces also ensure thermal equilib- 
rium even when dissipation is present, and in fact all 
our discussion about the quantum case is restricted to 
systems in thermal equilibrium. Of course we can also 
include a Langevin force in a classical model in order to 
balance dissipative effects and obtain an average energy 
consistent with thermal equilibrium. 



III. QED ENERGY DENSITY IN A UNIFORM, 
DISPERSIVE, NON-ABSORBING MEDIUM IN 
THERMAL EQUILIBRIUM 

One major distinction between the classical and QED 
theories, of course, is that in QED there is a nonvanish- 
ing zero-point energy associated among other things with 
Casimir effects. For purposes of comparison with results 
obtained in the following section when absorption is in- 
cluded, we reproduce here the thermal equilibrium QED 
energy density in a uniform, dispersive medium in which 
absorption is negligible. To do this we simply regard 
Eqs. (0) and ([8]) as expectation values when expressed in 
symmetrized form in terms of the quantized fields E and 
H. (We use circumflexes to designate operators.) For the 
zero-temperature (vacuum) state |0), for example, we use 
familiar expectation values, e.g., 

(0|E«(r,w) ■E ( - ) (r,w')|0) = |E(r, u)\ 2 5(uj - u'), 
<0|E ( -)(r,w)-EW(r,u/)|0) = 0, (35) 

where E^" 1 ") (r, w) and E^ - ) (r, lj) are respectively the pho- 
ton annihilation and creation parts of E(r, to), and obtain 
straightforwardly 
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by manipulations similar to those used in the preceding 
section. For the case of a uniform medium with negligible 
absorption, these expressions imply the (infinite) zero- 
point energy density 



ue + um 
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where A = 1,2 denotes polarization components. The 
zero-point squared amplitudes for the quantized field in 
a non-absorbing medium are pd| 
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(39) 



where nn(uj) = (efl(w)/Ltij(w)) 1 / 2 is the refractive index. 
Therefore we have the familiar result 



Ue + u M 



h f°° d 

duiuj 3 nj i (u>) — [wriR (id)] 



2tt 2 c 3 j 

h 

2^V 



duKV 3 'n 2 R (w) 



1 



V 9 (UJ) 



— - I dkk 2 uj 

2 f 1 
— / d^k-huj, 



8- 



where we have used the relation k = nn(uj)u)/c and the 
definition v g (oj) = c[d(nfjUi) / diu] -1 of the group velocity 
at frequency lu. Whereas we have obtained this result 
for the zero-point energy without taking absorption into 
account, we will show in the following section that it is 
valid in general for a uniform absorbing medium. 

The generalization to finite temperatures is similarly 
straightforward and yields, of course, 
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where the spectral energy density [ll|, 
n 2 R (uj)hw' 3 (I 1 
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IV. MODEL FOR A DISPERSIVE, ABSORBING 
DIELECTRIC MEDIUM IN THERMAL 
EQUILIBRIUM 

Following the work of many others [T3|, we model a 
dielectric medium as a collection of harmonic oscillators. 



Aside from the need to introduce oscillator strengths in 
order to obtain correct numerical results, the oscillator 
model is an excellent approximation if the atoms of a 
dielectric medium remain with high probability in their 
ground states. Each oscillator atom has a mass m and 
a natural frequency luq and is coupled to a reservoir of 
other harmonic oscillators responsible for the damping 
of its oscillations and homogeneous line broadening of 
its (electric-dipole) transition. The Hamiltonian for this 
model, including the electromagnetic field and its cou- 
pling to the atoms, is 



H = 



+ -muj 2 x 2 ^j + J dwfiw 



b H-b» + 



poo 

: / dujK.(w) y~}xj ■ [bj(u>) - b](w)]. 
Jo 



(43) 



The first two terms correspond in standard notation to 
the energy of the electromagnetic field, the atom oscilla- 
tors, and their coupling via the vector potential A(rj), 
where denotes the position of the jth atom. The third 
and fourth terms represent respectively the energy of the 
reservoir oscillators and their interaction with the atoms. 
The reservoir oscillators satisfy the bosonic commutation 
relations 

(40) [b ifi (uj), b] v {(j')] = SijS^Siui-J), b jv {J)] = , 

(44) 

where we use Greek letters to denote Cartesian compo- 
nents of vectors. The atom-reservoir coupling constant 
is chosen to be 



AM 



7T J 



1/2 



(45) 



in order that each atom's oscillations be damped at the 
rate 7, as shown below. 

From (HU and [ Xijj,,Pj V ] = ihSijSfiu we obtain in the 
dipolc approximation the Hciscnberg equations of motion 

x,-+wnx, = — E(r 7 -)H — / dwA(w)[b,-(w,i)-bt(w,i)l, 
m m Jo J 

(46) 



bj(uj,t) = -iujbj(Lj,t) + ^A(w)xj . (47) 



Using the formal solution of (|47l) in (|46l . it follows that 
ij+ufa = — E( ri ) + — F i -(f) + -L 

pQO pt 

x dujA 2 (uj) d£'xj(t')[2isinw(t'-t)], 
Jo Jo 

(48) 
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where the "Langcvin force" operator F^j (t) acting on the 
jth atom is 



dojA(uj) [bj (u, Q)e~ iut - b] (u, 0)e lLJt ] . 



(49) 
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The third term on the right-hand side of fl8l is 



mh 
27 



= 2 7 ^' dt%{t')-^8{t> -t) = - 7 %(t) . (50) 

We have omitted a divergent frequency shift which, when 
the atom-reservoir coupling is modified by a form factor 
to produce a finite expression, can be assumed to be in- 
cluded in the definition of the atom's transition frequency 
luq. Then (I48p simplifies to a "quantum Langevin equa- 
tion" [3]: 

... e - 1 » 

ij + 7% + w 2 Xj = — Efa) + —F Lj (t). (51) 
m m 

In the absence of coupling to the electromagnetic field 
we have, for times t 3> 7 -1 , 



(t) = / dwwA 2 (w) 



bj(u)e- 



LUq — UJ Z — IJUJ 



(52) 



Since we are working in the Heiscnbcrg picture, the ex- 
pectation value is over the initial state of the coupled 
system of oscillators. If wc assume that the reservoir is 
in an initial state of thermal equilibrium at temperature 
T, while the oscillator coupled to it is in its ground state, 
then 

= (hn(u)bj V (u')) -5 lJ 8 tnJ 8{uj -u/) 
1 



e hiu/k B T _ i l 3 u H 

A/ r (w)5y5 (tit/ 5(w - a/) 



(55) 



and 

z 1 ~ 2 1 2~2x 



3ft 7 f°° oj(uj 2 + oj 2 )jV(uj) 
duj 



3hj 



du> 



(oJq — UJ 2 ) 2 + J 2 LU 2 

oj(uj 2 + OJ 2 ) 

2\2 T 7a, ,2' 



(wg-w a ) a +7 



(56) 



The first term on the right is just the energy of the oscilla- 
tor in thermal equilibrium, and has a rather complicated 
closed form [13]; it becomes just 3hoj /[e huJo/kBT - 1] in 
the weak-coupling limit (7 — ► 0). The second term is the 
zero-point energy of the oscillator, that for u)\ > may 
be written as 



3H 

0J\ cos 

7T 



7 
2u ( 



3^7 / 

27T \ CJq 



(57) 



where ui c is a high-frequency cutoff [18j. It reduces to 
(3hujo)/2 in the weak-coupling limit. 



(We now write bj(w) in place of bj(w,0).) Similarly, 
using (j4"4"| . we obtain [l5[ 



2^7 



dww 2 cosw(i' — t) 



J (lj 2 - w 2 ) 2 + 7 2 w 2 



COSU>i(t' — t) 



2oj- 



■ sinwijf' — t\ 



-l\t'-t\/2 



(53) 



A. Noise Polarization 

The Heisenberg equations of motion for the electric 
and magnetic fields that follow from the Hamiltonian (|4"3"|) 
have, of course, the same form as their classical (Maxwell) 
counterparts 



V x E = 



15B 

'c at"' 



where us\ = [uj 2 — 7 2 /4] 1 / 2 ; thus the canonical commu- 
tation relation at equal times [xi^{t),Pj V {t)] = ihSijS^ 
is preserved in the coupling of the atom to the reservoir 

The energy expectation value of a single oscillator 
without coupling to the electromagnetic field is found 
similarly to be 



/I ~ 2 1 2-2\ ^7 I 



lo(lo 2 + UJ 2 ) 



II 



(oj 2 - UJ 2 ) 2 + J 2 0J 2 



x£[2$>)*i>)> + l]- 



„ TT 4tt» 1SE 
VxH=— J + - — , 

C COT 

that must be supplemented with 

V-B = , V-D = 0, 

where 

D = E + 4ttP, 

dP(r,t) 



J(r,i) = 



OT 



(54) 



P(r,t)=e^x J (i)5 3 (r-r J ), 



(58) 



(59) 



(60) 
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and, because our model does not induce any magnetic 
activity, we have B = H. It is advantageous to work in 
the frequency domain, so we write 

/■OO 

E(r,t) = / dw[E(r,w)e-^ t + E t (r,w)e^ t ], 
Jo 

/>oo 

H(r,t) = / ^[H(r,w)e- lwt + H t (r,w)e^ t ], 
Jo 

/■OO 

P(r,t) = / (L[P(r 1 w)e- <art + P t (r,w)e <<J '], (61) 
Jo 

where the Fourier transform of the polarization density 
may be written as 



P(r,w) = e^x.H^r-r,), 



(62) 



*,-(*) = / tL;[x J -(w)e- <w '+xt(a;)e iu *] . (63) 



It follows from l|5ip that 

P(r, W ) = e a /m )g3(r _ } 

^^e^m 3 



(jg — o;^ — 270; 

TVe 2 /™ », 
-5 2 • E(r,a;) 

0^5 - or - ijuj 

iNe/m ~ 
-2 5 : A c (w b r,w , 

075 — CJ^ — IJUJ 



(64) 



in the approximation in which we assume the atoms are 
continuously distributed with a density N and A c (to) = 
yj p m h"fuj/TT, with p m = m/N. 

From Maxwell's equations (|58l [59)) and the fact that 
V • b(r, w) = [lj|, we obtain 



V 2 E(r,w) + — E(r,w) = -4tt— P(r,w) 

c c 

-47riVe 2 /m w 2 -, 

-5 j = 2- E ( r ' w ) " 

— uj* — 170; c 

AirieN/m uj 2 - 
— 2 — : jAwbr,w, 

— UJ Z — IJUJ C Z 



(65) 



V 2 E(r,u;) + — e(,)E(r, W ) - __ K (r,o;), (66) 
where the complex permittivity is 



. , „ 4:TTNe 2 /m 
s(uj) = 1 - 2 2— — = 1 

= e R (uj) + iei(u). 



uj 2 — uJq + ijuj uj 2 — uj 2 + i^fuj 



(67) 



We have also defined the "noise polarization" at fre- 
quency uj: 



K(r,w) 



AmNe/m 
uj 2 uj 2 170; 



-A(w)b(r,w). (68) 



This contribution to the polarization arises from the 
Langevin force Fij(i) in the quantum Langevin equa- 
tion (|51[) . Its principal properties for our purposes are 
the thermal equilibrium expectation values 

(^(r,a;)) = (^t(r >W )) = > 
(i^(r,u;)it(r>')) = (kl(r,Lu)Ki(r\u')) = 0, (69) 

and 

(£t(r,o;)it(rV)) 



4fc / (w)(5^ ( 5(w - uj')8 3 (r - r') 
1 



(70) 



(i^(r,o;)£t(rV)) = 4Hei(u)S^5(u - uj')S 3 (t - r') 



1 



?hu)/k B T 2 



(71) 



all of which follow from ([55)1 and (0^)6^^(0-1')) = 0. 
Eqs. ([70l[7T|) constitute nothing else than the fluctuation- 
dissipation theorem, that we derived from the fundamen- 
tal assumptions of a canonical bath and a linear coupling 
to the matter. 

We can proceed formally now as in Lifshitz's paper Q 
and define operators g x (k, uj) by writing 



K(r,w) = d 3 k ^ <?A(k,w)e kA e 

" \ 1 '-> 



(72) 



A=1.2 



The solenoidal character of b(r, uj) implies directly in 
V • K(r, uj) = and therefore we can choose the vectors 
e k A such that k- e k A = 0, ekA • ekA' = 0, A = 1,2; we also 
take the ekA to be real. Then 



<?A(k,w) 



d 3 r K(r,uj) ■ e kx e 



-?kr 



= (^) / d 6 rKx{r,uj)e 



-ikr 



(73) 



and Eqs. ([55]) and (T4"4"l) imply the commutation relation 

MKoj),9{,{V,uj')] = ^e!(uj)6xx>6(uj-uj') 

x<5 3 (k-k'), (74) 

where again we make the uniform continuum approxima- 
tion for the spatial distribution of the material oscillators. 
Finally it will be convenient to introduce the operators 

Cx(k,u>) = [fe / (^)/2 7 r 3 ]- 1 / 2 3A(k, W ), (75) 

satisfying 

[Cx(Kuj),Cl,(k',uj')] =6xx>5(uj-uj')S 3 (k-k'). (76) 
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B. Remarks on the Lifshitz Theory 

It seems worthwhile as an aside to compare the formu- 
lation presented thus far with the Lifshitz theory. For 
this purpose we write 

\ (Kl(r, u)K„(r', J) + k„(T, oj)kt(r', u>')) 



Ahei{u)5^5(u - cu')S 3 (r - r') 



1 



(77) 



The right-hand side is equivalent to that in equation (1.2) 
of Lifshitz's paper Q , which in our notation has the form 



1 



1 



e hu/k B T _ i 2 



(78) 



the ((■■•)) denoting a classical ensemble average. This 
expression in the Lifshitz theory is a statement of a 
fluctuation-dissipation relation that Lifshitz attributes to 
Rytov 0. The difference between ([77]) and (JTSJ) reflects 
the fact that in the Lifshitz theory the thermal equi- 
librium electric, magnetic, and noise polarization fields 
are treated in effect as classical fluctuating fields; the 
fluctuation-dissipation relation (f78|) is used to relate the 
average of the square of the noise polarization to the 
imaginary part of the permittivity. The constant h ap- 
pears in ([75)) as a result of fixing the right-hand side 
such that the average over the classical ensemble for the 
squared fields matches the corresponding quantum ex- 
pectation values. In our (quantum) formulation based 
on the quantum Langcvin equation for the material os- 
cillators, the only nonvanishing contribution to the ex- 
pectation value of the square of the noise polarization at 
T = 0, for instance, is 

(K^r, u>)kt(r', w')> = 4fo 7 H<V^ - u/)<5 3 (r - r'), 

(79) 

twice the corresponding result in Lifshitz's paper. But 
Lifshitz's averages for the squared fields are the same 
as our corresponding quantum expectation values be- 
cause in his formulation both ((K*(r, u))K(r', a/))) and 
((K(r, u))K*(r', a/)}) contribute (equally) to these aver- 
ages. Thus the same Casimir force will be obtained in 
either approach because they both involve the same av- 
erage zero-point energy per mode, although of course the 
averages in the two approaches are fundamentally differ- 
ent. In the Lifshitz theory, in which forces between bodies 
are calculated using the stress tensor, there are no quan- 
tized fields, and averages of components of the stress ten- 
sor are over classical ensembles of stochastic fields, their 
statistical properties being determined by imposing the 
Rytov fluctuation-dissipation relation ([78[) . 

For the calculation of the Casimir force between per- 
fectly conducting plates, for example, a stochastic elec- 
trodynamical (SED) theory yields the correct force when 



h is introduced by requiring that there is a ze ro-p oint 
field energy (l/2)/iw per mode of frequency u) [2fl|. In 
SED, as in the Lifshitz theory both E* (r, ui) ■ E(r, ui) and 
E(r, ui) ■ E* (r, u>) contribute to the average of the squared 
electric field at T = 0, whereas in our quantized-field ap- 
proach only E(r,w) ■ W(r,ui) contributes. The Lifshitz 
approach to the calculation of Casimir forces may be re- 
garded as an application SED in which dissipation as well 
as finite thermal equilibrium temperatures are treated. 



C. Electric and Magnetic Fields 

An expression for the quantized electric field in an ab- 
sorptive dielectric now follows directly from Eqs. (|61jl . 
d66|), (|Z2J, and CfJ]): 

2 /„2 



E(r,i) = J d 3 k^2j^ duj^/he I (uj)/2ir ? - 



u> 2 /c 



k 2 — e(u))ui 2 1 'c 2 



CA(k, W )e kA e- i(wt - k - r >+h.c.. 



(80) 



From VxE = — (1/ c)<9B /dt we also obtain an expression 
for the magnetic field (reminding that /.t = 1 under our 
assumption in this section of a nonmagnetic medium): 



H(r,t) = ic f d 3 kY] f duj y / hei(uj)/2TT 3 
J x Jo 



k 2 — e(u))ui 2 1 c 2 
x lo' 1 (k x e kA ) e-^*- k - r) + h.c. (81) 

These expressions have the same form as the correspond- 
ing ones obtained by Huttner and Barnett [(| by Fano 
diagonalization of the entire system of coupled harmonic 
oscillators. 

One quantity of interest is the zero-temperature ex- 
pectation value of E 2 (r,£), for which the considerations 
above yield 



/>oo 

/ duiuj^eiiuj) 
c A Jo 



to ej(ui) 



e(w)w 2 /c 2 | 5 



4:Trk 2 dk 



/„ [k 2 -oj 2 € R {uj)/c 2 ] 2 +oj i € 2 {uj)/c i 
= doJuj 3 n R (uj), 

(82) 

where we have used the relations 6r(u>) = n 2 R {ui) — n 2 {uj) 
and ej(w) = 2n R (u))ni(uj) for the real and imaginary 
parts (n R and rij) of the refractive index. We note that 
this is the same form one would obtain by quantizing the 
field in a dispersive and non-absorbing medium, assuming 
a purely real permittivity e R (ui) [ll[. 
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For H 2 (r, t) we obtain the zero-point expectation value 



<H 2 (r,i)> 



Airk^dk 



o \k 2 ~t{io)uj 2 /c 2 \ 2 ' 



d 3 fc 



tu 2 e(uj)/c 2 



(83) 



The integral over fc diverges. However, it is obtained 
in the approximation that the atoms of the dielectric 
form a continuum, an approximation that is invalid when 
ka 1, where a is a typical interatomic spacing. To 
apply our macroscopic approach based on the charac- 
terization of the medium by a permittivity t{uj) — a long- 
photon-wavelength approximation implicit in the Lifshitz 
theory — we must "regularize" the integral (|83|) to extract 
a finite result. In this case it is convenient to introduce 
a Lorentzian cutoff to the fc-integral [2l[ and then, by 
taking advantage of the integrand even parity, a simple 
application of the residue theorem gives (in the contin- 
uum limit) 



d 3 k 



1 



P-2o J k 2 — euj 2 /c 2 1 -I- k 2 a 2 



2n 2 



+ 27 r 2 »-e 1 / 2 , (84) 



where, for clarity, we omitted the u-dependence of e(w). 
Finally, a direct substitution of (f8~4")l into ([83)1 gives 



^<H 2 (M)) 



E 



^+ W Rc6 3 / 2 
a 



(85) 



which, of course, diverges as a — > even before the tu- 
integration. We shall see in the next subsection that 
the divergent part is actually canceled out when all the 
contributions for the energy are taken into account. 



D. Energy Density 

To obtain an expression for the total energy density in 
the dielectric medium we start from Poynting's theorem 
in the conventional notation (symmetrized Poynting op- 
erator S = (c/8ir)[E xH-HxE]) and take expectation 
values over the initial state of the system consisting of 
the field, the dielectric atoms, and the reservoir: 







<9D 


dT> 


8ir] 


At" 1 


' ~dt 






<9H 


dH 


8^ J 


~dt ' 





'S) • xida 



According to the usual interpretation, the l.h.s. of ([86]) is 
the energy flux through a given surface S and, given that 



~E)dV 



H)dV. (86) 



we are assuming thermal equilibrium within our system, 
it should vanish. Thermal equilibrium also allows us to 
identify the rate of change with time of the expectation 
value of the total energy W per unit volume: 



dW 



1 

8tt' 



E 



<9D 

~dt 



<9D 

~dt 



E) 



J_d_ 

8^dt 



H 



(87) 



which, of course, also vanishes. For the system under 
consideration D = E + 47rP e + K, where P c is the part 
of the polarization giving rise to the dielectric permittiv- 
ity e(u>) and K is the noise polarization defined by 
Thus D = D £ + K and 



dW 
~dt 



j i 

dt 



dt 



where we define 
dWi 



dt 



and 



— (E 



dW 2 



dT> f dT> e 



dt 



dt 



E) 



1 a m 2 



dt 



i 

8^ 



E 



dK 

~dt 



dK 
~dt 



E) 



(89) 



(90) 



Before proceeding with the calculation of W we note 
the following identity that follows from our model of the 
dielectric: 



dW 
~dt 



'dt ^ 



1 A2 

-mx j 



<5 3 (r-r,) 



E + H z 



j 

J_d_ 

I^di 



Yfey^rr&j) -%-F ii ]5 3 (r-r J -)). 



(91) 



The first term is the rate of change of the energy density 
(kinetic plus potential) of the oscillators constituting the 
dielectric, and the second term is the rate of change of 
the energy density of the electromagnetic field. In the ab- 
sence of any dissipation (7 = and therefore Flj = 0), 
the third term on the right vanishes, and W = u is just 
the total (matter-plus-field) energy density. The third 
term accounts for the effect of tjhe reservoir on the di- 
electric oscillators: 27^ J (|mx 3 )5 3 (r — r^) is the rate 
of change of kinetic energy density due to the dissipa- 
tive effect of the reservoir, while ^ ■ Xj • FLj8 3 (r — r 3 ) is 
the rate of work per unit volume done by the Langevin 
forces on the dielectric oscillators. In the absence of the 
electromagnetic interaction these effects cancel, and the 
third term in (f9Tj) again vanishes. There is a close formal 
similarity between (|9T|) and the corrresponding expres- 
sion (|34[) that follows from the classical oscillator model. 
The essential physical difference between (|9Tj) and (|34"|) 
lies simply in the effect of the Langevin force term in the 
quantum-clcctrodynamical expression of energy conser- 
vation. 
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To obtain the total energy density we focus first on 
the case of zero temperature, as the result for finite tem- 
perature requires only a simple extension of the zero- 
temperature calculation, as discussed below. Using (j6"Tj) 
plus 



<9D e 
~dt 



dww[e(w)E(r, uo)e 



- £ *( W )Et(r )W )e +farf ], 
and integrating over t, we obtain 

lo' 'e*(o/) — loe(lo) 



(92) 



1 /■CO /■OO 



x (E A (r,u)-E),(r,a.'))e- , ("-"')' 



where we have used the fact that the vacuum (zero- 
temperature) expectation value (E A (r, w)-Ea'(i", lo')) = 

while (E^(r, w) • Et,(r, u>')} vanishes unless A = A' and 
lo = lo' . To deal with what appears to be a singularity at 
lo = lo' we rewrite (|93[) as a sum of two identical terms 
and interchange lo and lo' in the second one, to get 



;H 2 (r,t)), 



(93) 



8t jio Jo W - lu 

1 ST fZ ' fZ < ' I 'u f " (EA(r, w) ■ E* (r, q/)) e -i("-"')« - (E A (r, a/) ■ E* (r, a^e^^ 

+ ^(H 2 (r,t)). (94) 



Next we use (|80|) and (|T6[) to write the vacuum expectation value 



(E A (r,^).Et(r,^)) = (EA(r, W ')-Et(r, W )) = A e/M ^ j d ^ k __l___ S (co - <J), (95) 



that allows us to readily evaluate the first term in (|M|) by noticing that 

u/e R (u/) - loc r (lo) d 

hm = — we fl w . (96) 

w'^w w' — w aw 

The second term is calculated by realizing that the zeroth order contributions in (lo — lo') in the numerator cancel 
each other, while the first order terms produce a contribution linear in the elapsed time t 

p — i(uj—uj )t p i(uj— uj )t 

lim —— r = it. (97) 

ul'^u) 2{L0 — LO) 



Therefore 



w ^ i > = 5,gi:r*(s M+ ' ti, ')' A '/ rt if-^f + > M 

/■<» d 1 fi, Z" 00 

E y o fa3n ^ N + to (fi2(r '' )) + i, 4^E y o ^ 4 ^/ (98) 



8tt 2 c 3 



r 



where we have again used the relations €r = n R — nj change in time of W\ (t) is a positive constant, given by 
and ej = 2n R nj and we are leaving the w-dependence the last term in (|98p l22j|. This implies that this term 
implicit in both e and n. Let us note that the rate of is responsible for creating heat in any given volume V 
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(meaning an inward flux of energy). Since we know that and use (|72j) . ([75]) . and ((80)) to relate K,\(k, uj) and 
in thermal equilibrium the total flux should vanish, this E A (k, uj): 
energy increase must be balanced out by an energy de- 
crease coming from W-z(t). 

To evaluate W2 as given in Eq. (|90|) we first define 
K(k, uj) by writing 



K(r,£) = dtu J ti 3 £;VjK A (k, 



uj)e~ tuJt ^ r 



Kt(k, W )e^e- ik - r ], 



K A (k, uj) = ^[k 2 -e(uj)uj 2 /c 2 }E x (k,uj). (100) 



Then, after inserting ([99]) and (fTOOf in (J90J) and a few 
(99) algebraic steps, we get 



W 2 (v,t) 



16tt 4 c 2 



Jo 



J2 / dw' / - v /e / (w)£j(w')*(w-w') / d 3 fc 



a i{uj-uj')t 



k 2 — e(uj)uj 2 I c 2 k 2 — e* (uj)uj 2 / c 2 

(101) 



and, proceeding as in the evaluation of Wi, we obtain 

/■« ,.,2,. 

^E Rc , 
/0 



W 2 (r,t) - ^VRc / duj' / <L; ye/Me/Cw'Ww-u/) / d 3 fc 



k 2 - e(uj')uj' 2 /c 



t 



Att 2 c 3 



poo 
x JO 



(102) 



r 



where we have used the integral / d 3 k\k 2 — euj 2 /c 2 \~ 2 = 
2TT 2 cn R / ejuo (as in Eq. (|82p) to obtain the second term. 
Now we see clearly that the time dependent term in W2 (t) 
precisely cancels the one in Wi(t), ensuring thermal equi- 
librium. Let us note also that the first term in (|102p con- 
tains the same fc-integral as the one present in (|83p . and 
therefore we may invoke Eq. (|84|) to evaluate it. The 
first term of (|102|) is then equal to 
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8tt 4 c 2 ^ 



duj' I dui 



uj 2 uj' 



2tt 2 



S(uj-uj'), 
(103) 

where we again left implicit the uj- and o/-dependences in 
e/. The apparent singularity in the uj' — > uj limit may be 
dealt with by using the procedure described just before 
Eq. (|94|) . and then after some trivial steps we get 



8ir 2 c- 



E 



dujuj 2 — , 



(104) 



what cancels exactly the first term in (|85[) . We still have 
to work on the contribution of the second term of (1841) to 



the first term of (|102p . that leads to 

h 



duj I duj' 




x lim 

Uj' — >UJ 

h 



UJ — UJ 

r°° d 

3 Vim / dujuj 2 e 1 (uj) — [uj 2 e 1/2 (uj)}. (105) 
, Jo dui 



The total energy density is obtained by adding (|98|) and 
(fT02jl and using (fT05l) . 185]) : 



8tt 2 c 3 



E 

A 



dww 3 ^ Re 



n^ — (we) + e 3 / 2 
duj 



uj duj 



(106) 



The above expression of the energy density of a uniform, 
dispersive and absorbing medium, is the most important 
result of this paper. Using e(uj) = n (uj) and the follow- 
ing relations 

nR ~ d ~ J ( u}eR ) = ( n « ~ n2 i) n R 

(r> dTl R o dUl 

+ UJn R 2n R — 2nj — — 

\ dw duj 

Ree 3/2 = (n 2 R - n 2 j)n R - 2n R n 2 , 

ej d , 2 n . 2 r, dni 

Im— (uj v e ) = 4:n R n r + 2n R njuj— — 

duj 



uj duj 



(107) 
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and summing over polarizations we obtain our final ex- 
pression for the vacuum expectation value of the total 
energy density: 

which is just 140(1 : the QED zero-point energy density 
depends in exactly the same way on the refractive index, 
regardless of whether absorption is accounted for, and 
in fact it depends only on the real part of the refractive 
index. 

The same conclusion holds for finite temperatures. 
In this case both (E^(r,w) ■ Ey(r,u}')) and (E A (r,w) ■ 
E A ,(r, u/)) make nonvanishing contributions to W\ and 
W 2 . Using 

(Cl(k, u)C x , (k', J) + C A (k, u)C{, (k', «')) 

= 5x, x .5{w - cu')S 3 (k - k') coth Jf±=, (109) 

we arrive straightforwardly at exactly the formula (pfTj) 
for the total energy density in thermal equilibrium. 

V. CONCLUDING REMARKS 

We have shown that the approximate classical expres- 
sion (TT]) for the energy density in a band of frequencies at 
which absorption in a dielectric can be ignored is in fact 
exactly correct as a spectral average value in (i) classical 
theory in the case where the fields at different frequen- 
cies within the band undergo uncorrelatcd fluctuations, 
and (ii) QED at zero temperature or more generally at 
thermal equilibirum. 

Using the model of a dielectric medium as a collec- 
tion of harmonic oscillators, and including the coupling 
of each oscillator to a reservoir of oscillators that give 
rise to dissipation and a Langevin force on each oscilla- 
tor, we have shown how a noise polarization results from 
these reservoirs and compared it with that employed in 
Lifshitz's well-known theory of Casimir effects. From 
this comparison we concluded that the Lifshitz theory 
is actually a classical stochastic electrodynamical theory. 
We arrived at quantized electric and magnetic fields hav- 
ing the same form as in the Huttner-Barnctt approach 
in which the complete Hamiltonian is diagonalized, and 
showed that the expectation value of the total energy 
of the system of dielectric oscillators, reservoirs, and the 
electromagnetic field has the same form in thermal equi- 
librium, including the limiting case of zero temperature, 
independent of whether we take dissipation into account 
in quantizing the field. 

Our treatment allowed us to derive the celebrated 
fluctuation-dissipation theorem, and also to show explic- 
itly that it ensures that in thermal equilibrium the total 



energy of the system of oscillators, reservoir and electro- 
magnetic field is constant in time. When absorption is 
present, there is a coupling between the system of oscil- 
lators and the reservoir and an energy exchange between 
them. In our example of dielectric medium modeled by 
a collection of harmonic oscillators, we have explicitly 
shown that a positive energy rate W\ > arising from 
the interaction of the EM field with the system is exactly 
canceled by a corresponding negative energy rate com- 
ing from the interaction of the system with the reservoir, 
W<i = —W\ < 0. We should stress that this energy rate 
balance is absolutely general, and applies not only to 
absorbing dielectric media, as treated here, but also to 
arbitrary dissipative materials, e.g. metals modeled by a 
dissipativc Drudc permittivity. This is merely a manifes- 
tation of the fluctuation-dissipation theorem, and as such 
holds for Lifshitz theory, which assumes thermal equilib- 
rium. 

These conclusions cannot be applied directly to the 
question raised in the Introduction: why do derivations 
of Casimir forces that start from calculations of changes 
in zero-point energy of presumed dissipationless media 
appear in the end to produce the same results as in the 
Lifshitz theory, where a correlation function of a fluctu- 
ating ( "noise" ) polarization is related to the imaginary 
(dissipative) part of the permittivity? To address this 
question we must take into account a hallmark of Casimir 
effects, namely the role of boundaries. We have restricted 
ourselves here to the model of perfectly uniform media 
with no boundaries; the role of boundaries will be ad- 
dressed in a forthcoming paper Q ■ 

Finally we mention two examples, not involving bound- 
aries, where these conclusions are consistent with known 
results. The first example is spontaneous emission of 
an atom embedded in a dielectric medium. The Ein- 
stein A coefficient for the rate of spontaneous emission at 
the electric dipole transition frequency ujq is proportional 
to the zero-temperature expectation value of E 2 (r, ojq), 
which from Q82|) is seen to be proportional to ur^uo). 
Therefore the A coefficient for an atom in a dielectric 
medium with complex refractive index nji{uj) + inj(tu) is 
just tir(ujo) times the free-space A coefficient [23j | . This 
assumes the continuum approximation for the dielectric; 
near-field interactions of the embedded excited atom with 
host atoms, including local field corrections, result in a 
rate of energy loss by the excited atom that depend on 
nj(wo) 0. 

The second example, less straightforard but more 
closely related to Casimir forces, is the van der Waals 
interaction between two atoms embedded in a dielectric 
medium. It has been shown that the interactions between 
electrically or magnetically polarizable atoms can be ob- 
tained from the quantized electric and magnetic fields 
in a nonabsorbing medium [25} . Absorption affects the 
final expression for the interaction only after the permit- 
tivity (or permeability) is properly regarded as a complex 
function of frequency, analytic in the upper half of the 
complex frequency plane. As in the Lifshitz theory, and 
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as in the present work, the calculations leading to this 
result are based on the continuum model of the dielectric 
medium. 
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